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1. Introduction 



The two-level pairing model describes a finite system which undergoes a second- 
order quantum phase transition between weak-coupling and strong-coupling dynamical 
symmetry limits. This quantum phase transition is characterized by singularities in 
^ ! the evolution of various ground state properties as the pairing interaction strength is 
i varied: (1) a discontinuity in the second derivative of the ground state eigenvalue, 
(2) a discontinuity in the first derivative of the quantum order parameter, which is 
defined by the relative population (A''2) — (Ni) of the two levels and is analogous to the 
magnetization parameter in the Ising model, and (3) a vanishing energy gap A between 
the ground state and first excited state with the same conserved quantum numbers, and 
thus a singular level density p ~ A~^. Although true singularities in these quantities 
only occur in the limit of infinite particle number, "precursors" are found at finite N, 
which approach the singular limit according to definite power- law scalings [1-8]. The 
quantum phase transition in the two-level pairing model has long been of interest for 
applications to nuclei [9-11]. It has recently served as a testbed for considering phase 
transitional phenomena, including the finite-size scaling just described, excited state 
quantum phase transitions [8,12-15], thermodynamic properties [16], decoherence [17], 
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and quasidynamical symmetry [18], as well as for developing the theoretical methods for 
treatment of these phenomena, including continuous unitary transformation [5, 7] and 
Holstein-PrimakofT boson expansion [19]. 

Finite systems with pairing-type interactions, and consisting either of bosons or 
fermions, occur in a broad variety of physical contexts. Fermionic examples include 
superconducting grains (electrons) [20] and the atomic nucleus (nucleons) [21]. Bosonic 
examples include the s-wave and d-wave nucleon pairs of the interacting boson model 
(IBM) [22], which themselves undergo a bosonic pairing interaction in the description 
of nuclear quadrupole collectivity, and condensates of trapped bosonic atoms [23,24]. 

The Lie algebraic properties of the two-level bosonic and fermionic systems are 
closely parallel, but the differences which do arise fundamentally affect the irreducible 
representations (irreps) under which the eigenstates transform and therefore are essential 
to defining the spectroscopy of the system. Two complementary algebraic formulations 
are relevant to the description of finite pairing systems [25-30]: a unitary algebra is 
spanned by the bilinear products of a creation and anniliation operator [31,32], and a 
quasispin algebra [25,33,34] is defined in terms of creation and annihilation operators 
for time-reversed pairs of particles. These structures are intertwined by duality relations, 
in particular, relating irreps of the quasispin algebra with those of an orthogonal (in the 
case of bosons) or symplectic (in the case of fermions) subalgebra of the full unitary 
algebra. Such relations have often been used [30,35-40] to effect simplifications of the 
calculations for two-level and multi-level systems. 

In this article, the duality relations between the unitary and quasispin algebraic 
structures for the two-level system are systematically established. In particular, 
attempts to compare results across two-level systems with different level degeneracies 
or between the bosonic and fermionic cases (see Ref. [14]) raise the question as to which 
differences in spectroscopic results are superficial, i.e., originating from an imperfect 
choice of correspondence between the Hamiltonian parameters for the two cases, and 
which are due to more fundamental or irreconciliable distinctions. Therefore, a main 
intent of the present work is to resolve the relationships between the disparate forms 
of the Hamiltonian which arise in the definitions of the dynamical symmetries and in 
numerical studies of the transition between them. These Hamiltonians include: (1) the 
Casimir form defined in terms of the unitary algebra, (2) the pairing form used in studies 
of the fermionic system, which is essentially defined in terms of quasispin operators, and 
(3) the "multipole" form traditionally considered for physical reasons in bosonic studies. 
The relations are established in a fully general fashion which uniformly accomodates 
arbitary level degeneracies (ni and 712), for both the bosonic and fermionic cases. 

Previous work on two-level systems has concentrated either on the so-called s- 
b boson models, with level degeneracies Ui = 1 and n2 > 1 {i.e., for which one 
of the levels is a singlet), or on fermionic models of equal degeneracies (ni = n2). 
The observations outlined here are intended to provide a foundation for more detailed 
future work, allowing for the most general choice of level degeneracies. The results 
are provided as a basis for algebraic studies of the quantum phase transitions, excited 
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state spectroscopic structure, and classical geometry [41,42] of two- level and multi-level 
pairing models. Although the discussion is presented for two-level systems, for the 
sake of clarity, many of the results carry over to multi-level systems essentially without 
modification. 

After a brief summary of the dual algebraic structures for the many-body 
problem in general (Sec. 2), the unitary algebraic structure is presented in detail, 
including categorization of the subalgebra structure, classification of the irreducible 
representations, construction of the generators, and identification of the Casimir 
operators, all in a unified form for bosonic and fermionic cases (Sec. 3). The simpler 
quasispin structure is also reviewed (Sec. 4). Duality relations are then established 
between the unitary (or Casimir) and quasispin (or pairing) formulations of the 
Hamiltonian (Sec. 5). These are explicitly related to the spectral properties of the 
two-level system through numerical calculations across the quantum phase transition, 
illustrating basic distinctions between the bosonic and fermionic cases, when calculated 
for bosonic and fermionic systems with with similar level degeneracies and/or similar 
particle number (Sec. 6). 

2. Bosonic and fermionic algebras 

The fundamental Lie algebra describing transformations of a many-boson or many- 
fcrmion system is spanned by the bilinear products of creation and/or annihilation 
operators al^/a]^, al^iCim, and a^/a^ {e-9-, Refs. [29,30]). For bosons, the resulting 
algebra is Sp(2n, M), and for fermions it is S0(2n), where m and m' — 1, . . ., n range 
over the single particle states of the system. Two important sets of subalgebras arise: 
number- conserving subalgebras and mimher-nonconserving (quasispin) subalgebras. 

The restriction to number-conserving operators, spanned by the elementary one- 
body operators al^,am, constitutes a U(n) algebra. The U(n) algebra contains a 
subalgebra SO(n) for the bosonic case or Sp(n) for the fermionic case. If each single- 
particle creation operator a]^ is associated with a time-reversed partner a^^, these SO(n) 
or Sp(n) subalgebras arc defined by the property that they leave invariant the "scalar" 
pair state X^^OmajnlO) [30]. This special property underlies the duahty relations with 
the quasispin pair algebra considered in the present work. More specifically, we consider 
rotationally-invariant problems, for which the single-particle states may be identified 
as the 2j + 1 substates of single-particle levels of various angular momenta j {i.e., 
j-shells, in the nomenclature of nuclear physics, which we adopt for either bosonic 
or fermionic levels). Then the creation operators are of the form aj„ — )■ a|,^ and 
at — )■ (—)""" ^""afc _m f*^^ level. For such rotationally-invariant systems, the SO(n) 

or Sp(n) subalgebras in turn contain the physical S0(3) ~ SU(2) angular momentum 
algebra. Although we follow the convention of denoting the angular momentum algebra 
S0(3) in the bosonic case and SU(2) in the fermionic case, there is no material distinction 
between the algebras. In general, there may also be other, intervening subalgebras in 
the chain. 
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Alternatively, the scalar pair creation operator S+ — IX^aj^a^, the scalar pair 
annihilation operator = ^^a^flm- and the number- conserving operator jS'o = 
iS(^m^'n + Oamdln) closc Under commutation, where ^ = + for bosonic systems or 
6 = — for fermionic systems. These operators define a number-nonconserving pair 
quasispin algebra, either SV({1, 1) for bosons or SU{2) for fermions. The caligraphic 
notation for the quasispin algebras is adopted [30] to avoid ambiguity between the 
SU{2) quasispin algebra and the SU(2) angular momentum algebra. 

In summary, for a bosonic system, the subalgebras under consideration are 

^P^^'^'^^^l 5W(1,1)DU(1), 
and, for a fermionic system, they arc 

^^^^""^ \ SU{2) D U(l). 

The subalgebras of U(n) are useful in the classification of states not only for the pairing 
Hamiltonian (defined in Sec. 5.2) but also for a much richer range of Hamiltonians [43]. 

A close relation between unitary chain and quasispin subalgebras arises since 
the quasispin and orthogonal or symplectic algebras may be embedded as mutually 
commuting "dual" algebras within the larger Sp(2n, M) or S0(2n) algebra. The algebraic 
foundations are discussed in detail in Refs. [25-30]. Here we simply note that duality 
denotes the situation in which the states within a space may be classified simultaneously 
in terms of two mutually commuting groups (or algebras) Gi and G2, such that the irrep 
labels of the two groups are in one-to-one correspondence. For the present problem, the 
embedding and associated labels are given, for the bosonic case, by 

Sp(2n, R) D [SO(n)D • • • dS0(3)] (g) [SU{1, 1)dU(1)] (3) 

[v] J S N 

and, for the fermionic case, by 

S0(2n) D [Sp(n)D • • • dSU(2)] (g) [5W(2)dU(1)]. (4) 

{v} J S N 

The seniority label v (Sec. 3.2) and quasispin label S (Sec. 4) are in one-to-one 
correspondence, i.e., specifying the value of one uniquely determines the value of the 
other, and vice versa. 

The duality relations hold equally well regardless of whether the n-dimensional 
single-particle space is construed to consist of a single j-shell {n — 2j + 1, odd for 
bosons or even for fermions), two j-shells {n — rii + 712), or, indeed, multiple j-shells. 
The one- level case has been considered in detail {e.g., Ref. [29]). However, we find 
that the detailed construction of operators for two-level systems within the context of 
this duality, as needed for spectroscopic studies of these systems, requires elaboration. 
Although, for simplicity, we consider only the case of two levels, the results may readily 
be generalized to additional levels. 
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3. Uniteiry algebra 

3.1. Subalgebra chains 

Consider the U(n) subalgebra chains for the two-level system, consisting of either bosonic 
or fermionic levels, of possibly unequal degeneracies. If the levels are j-shells of angular 
momenta ji and j2, the level degeneracies are ni = 2ji + 1 and n2 = 2j2 + 1, and the 
total degeneracy of the system is n = ni + n2- For the bosonic case, we have 

( S0(ni+n2) 
[jv] Ui(ni) (g) Uslns) 

I [A^i] m 

DSOi(ni) ® S02(n2)DSOi(3) ® S02(3)dSOi2(3) (5) 

[vi] [Vi] Ji J2 J 

and, for the fermionic case, we have 



U(ni + 712)3 



Sp(?^l + n2) 

M 

"{AT} '"^^^ ] Ui(ni) (g) U2(n2) 

{Ni} {N2} 

DSpi(ni) (8) Sp2(n2)DSUi(2) ® SU2(2)dSUi2(2), (6) 

{vi} {V2} -Jl J2 J 

where Ui = 2ji + 1 and n2 = 2j2 + 1- The irreducible representation labels, indicated 
beneath the symbol for each algebra, are defined in Sec. 3.2, and the algebras themselves 
are constructed explicitly in Sec. 3.3. Throughout the following discussion, bosonic and 
fermionic cases will be considered in parallel. 

The subalgebras summarized in (5) and (6) are generically present, regardless of the 
level degeneracies rii and 7^2, for rii and 7^2 > 2. However, several clarifying comments 
are in order: 

(1) The important special case of a singlet bosonic level {jk — 0) leads to rife = 1, 
and the corresponding orthogonal algebra SOfc(nfc) is undefined. The label may still 
be defined, in a limited sense, through the quasispin, as noted in Sec. 4. Two-level 
bosonic problems in which ji = are termed s-b boson models. These include the 
Schwinger boson realization (ji = j2 = 0) of the Lipkin model [44]. The subalgebra 
chains and labeling schemes for the s-b models were considered in Ref. [14]. 

(2) For a fermionic level with jk = ^, and therefore = 2, the symplectic algebra 
Spfc(nfe) in (6) is identical to the SUfe(2) angular momentum algebra. 

(3) Additional subalgebras of \J{ni-\-n2) may also arise, parallel to chains indicated 
above and still containing the angular momentum algebra, e.g., for the interacting boson 
model (ni = 1 and 712 — 5), there is a physically relevant chain U(6) D SU(3) D 
SOi2(3) [22]. However, since these chains are not directly relevant to the pairing 
problem and cannot be treated in a uniform fashion for arbitrary ni and n2, they are 
not considered further here. 

(4) Further subalgebras may also intervene between SO(n) and SO (3), or between 
Sp(n) and SU(2), the classic example being the appearance of the exceptional algebra 
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G2 in the chain S0(7) D G2 D S0(3) [31]. 

(5) Whenever two reahzations of the same algebra commute with each other, the 
sum generators Gi = G-^^ + Cf'' also form a new realization of the algebra, as in 
ordinary angular momentum addition. For instance, addition of the angular momentum 
generators for the two levels [L = L^^^ + L^'^^] gives the generators of the sum angular 
momentum algebra SOi2(3) (bosonic) or SUi2(2) (fermionic) in (5) and (6). However, 
more generally, if the two levels have equal degeneracies {rii — 712 = ^112), such a 
combination of generators may also be made higher in the subalgebra chains (5) and (6), 
yielding Ui(ni2)(8) 112(7112) D Ui2(ni2) D SOi2(3) and S0i(ni2) (8) 802(7112) D 8012(^12) D 
SOi2(3) for the bosonic case, or similarly Ui(ni2) ® U2(ni2) D Ui2(ni2) ^ 8Ui2(2) and 
8pi(7ii2) <S> 8p2(ni2) D 8pi2(ni2) D SUi2(2) for the fermionic case. 

3.2. Branching 

The branching rules for the irreps arising in the bosonic or fermionic realizations of the 
algebras in (5) and (6) provide the classification of states for the two-level pairing model. 
Some, but not all, of these branchings can be expressed in closed form. 

For the bosonic realization of U(n), the symmetric irreps [A^] = [A^O . . . 0] (with 
n labels) are obtained, where is the occupation number. For SO(n), the irreps are 
[v] = [vO...O] (with [n/2j labels). 

The U(n) SO(n) branching is of the type considered by Hammermesh [45] and 
is given by 

v = {Nmod2),...,N -2,N. (7) 

This rule applies both to the branching U(ni +n2) — ?■ SO(ni +712) and to the branching 
associated with each of the two levels in Ui{ni) 1)2(712) S0i(7Zi) ® 802(772). Note 
that n is odd for a single bosonic level and is even for the two-level system. [The 
U(7i) — > Ui(7ii) (8) 1^2(712) branching rule follows trivially from additivity of the number 
operators, N ^ Ni + N2.] 

For the branching 80(7ii + 712) — >■ 80i(7ii) 802(712), the allowed vi and V2 are 
obtained by considering all partitions of v as 

v^Vi + V2 + 2n^ (7i^ = 0,l,...,[7;/2j). (8) 

This rule may be verified by dimension counting arguments, that is, dim[t'] = 
^^^^2 dim[7;i] dim[7;2] , using the Weyl dimension formula [45] . Notice that for the bosonic 
system (in contrast to the fermionic case below) the branching rule for 80 (711 + 712) — > 
80i(7ii) <S> 802(712) is independent of the level degeneracies ui and 712, and the total 
occupation number N infiuences the allowed 80i(7ii) and 802(712) irreps only through 
the constraint (7) on v. 

The allowed partitions v — )■ {V1V2) are listed, for low v, in Table 1. As a concrete 
example, for N = 2, the allowed SO(ni + n2) irreps have v = and 2, with branchings 
to SOi(ni) (X> 302(712) given by the corresponding rows of Table 1. As a specific example 
of the equivalence of dimensions, consider the case of the two-level bosonic system with 
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Table 1. Partitions v {vi,V2) of the form permitted by the SO(ni +112) 
SOi(ni) (g) 802(^2) branching rule (8). The partitions permitted by the Sp(ni +112) — >■ 
Spi(ni) (g) Sp2(n2) branching rule (10) are a subset of these. 
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(2,0), (1,1) 
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(0,0) 







= j2 = 1, thus described by SO (6) D S0i(3) S02(3). The -y = 2 irrep of S0(6) has 
dimension 20, while the corresponding SO(3)(8)SO(3) irreps likewise have total dimension 
dim(2, 0) + dim(l, 1) + dim(0, 2) + dim(0, 0) = (5)(1) + (3)(3) + (1)(5) + (1)(1) = 20. 

The branchings of the form SO(n) — )■ S0(3), needed for SOi(ni) ® 802(712) — >■ 
S0i(3) ® 302(3), are more complicated and, in general, involve missing labels. However, 
such branchings occur widely in physical applications, and general methods exist for the 
solution based on weights or character theory [29,46]. An explicit multiplicity formula 
is obtained in Ref . [47] , applicable to the symmetric irreps arising in the present bosonic 
case (5). Finally, the reduction S0i(3) ^ S02(3) — )■ SOi2(3) is governed by the usual 
triangle inequality for angular momentum addition. 

The branching rules for the fermionic case are nearly identical, with a few 
modifications. For U(n) we obtain the a?^^^symmetric irreps {A^} = [1, . . . , 1, 0, . . . , 0], 
that is, with unit entries (out of n labels total). Similarly, for Sp(n) we have 
{v} = [1, . . . , 1, 0, . . . , 0], with V unit entries (out of n/2 labels total). 

The branching rule for U(n) — )■ Sp(n) [45] requires the modification of (7) to 

V ^ {N' mod2),...,N' -2,N', (9) 

where A"' = min(A", n — A"). Notice, therefore, that v < \n. This rule apphes to 
U(ni + n2) Sp(ni + n2) and to Ui(ni) ® 112(722) -> Spi(ni) ® Sp2(n2). 

The values of vi and V2 arising in the branching Sp(ni + 77,2) — ?> Sp]^(ni) (g) Sp2(n2) 
are again given by the partitioning condition (8), but now subject to an additional 
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Table 2. Branching Sp(4) — >■ Spi(2) <S> Sp2(2) for irreps of the two-level fermionic 

system consisting of two j = 1/2 levels. Allowed irreps are a subset of those listed 
in Table 1. The Cartan labels of these irreps with respect to the isomorphic chain 
S0(5) D S0(4) are also indicated. 



N V [Ai,A2]sO(5) {vi,V2) [Al,A2]sO(4) 









[0,0] 


(0,0) 
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2 


[1,0] 


(1,1), (0,0) 


[1,0], [0,0] 


3 


1 
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[0,0] 


(0,0) 


[0,0] 



Table 3. Branching Sp(20) — )■ Spi(lO) (g) Sp2(10) for irreps of the two-level fermionic 
system consisting of two j = 9/2 levels, with N = 10. The fermionic branching rule (10) 

restricts the allowxMl irrei)s. relati\(> to tlK)S(> listed in Table 1. for r > 4" 12 = 5. 



V ni2 
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(0,0) 
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(2,0), (1,1), (0,2), (0,0) 
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6 


4 
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2 


(5,3), (4,4), (3,5), (4,2), (3,3), (2,4), (3,1), (2,2), (1,3), 






(2,0), (1,1), (0,2), (0,0) 


10 





(5,5), (4,4), (3,3), (2,2), (1,1), (0,0) 



constraint, so 

v = vi+V2 + 2n^ {n^ = 0,1, . . . ,[v/2\) 
\{vi - V2) - l{ni - 712)1 < l{ni + 712) - V, 

as can again be verified by dimensional counting. The branching rules (9) and (10) 
together automatically enforce < |niandi'2 < 1^2 ■ For two levels of eguaMegeneracy 
{rii — 712 = "^12), the constraint simplifies to 

\vi — V2\ <ni2 — v (11) 

and only serves to exclude {vi,V2) values when v > \ni2- 

For illustration, branchings for the low-dimensional case Sp(4) — > Sp(2)(8)Sp(2) (two 
3 — \ levels) are given in Table 2. The chain Sp(4) D Sp(2) ® Sp(2) is isomorphic to 
the canonical chain SO (5) D SO (4) of orthogonal algebras, and the branchings given in 
Table 2 therefore also follow from the SO(n) canonical branching rule [48-50]. The SO (5) 
Cartan labels are given by [Ai, A2] = + ^'2)1 "^2)], where [A'^, A2] are the Sp(4) 
Cartan labels, and the S0(4) Cartan labels are given by [Ai, A2] = [1(^1 + ^2), \{yi — V2)] 
(see Ref. [51] for a summary of notation for this chain). Branchings for the higher- 
dimensional case Sp(20) D Sp(lO) ® Sp(lO) are given in Table 3. 
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3.3. Generators 

Of the generators for the algebras in chains (5) and (6), those involving a single level 
are well known [32]. Here we must construct the generators for the two- level system. 
Since the subalgebra chains terminate in the two-level angular momentum algebra, 
SOi2(3) ~ SUi2(2), it is most natural to express the generators as spherical tensors 
with respect to this angular momentum algebra. 

First, let us briefly review the results for the single j-shell, with creation operators 
a]^ and annihilation operators {'m = —j, —j + 1, • • •, j)- In the case of a single 
bosonic level, with degeneracy n = 2j + 1 (j integer), the subalgebra chain for U(n) 
is U(n) D SO(n) D • • • D S0(3) [see (1)]. The generators of U(n), in spherical tensor 
form, are the bilinears 

= (atx5)(^) ((7 = 0,l,...,2j), (12) 

where 7 = —g, —g + 1, • • ., +g. The product of two spherical tensor operators is dcflncd 
by (A" X B^y:^ = Ea/3(aa&/3|c7)A^5^. We follow the time reversal phase convention 
A'^ = (— )"~"yl"^ [52]. Thus, e.g., = {—y~"^a_rn, where the time reversal phase 
factor is required for the annihilation operator to transform as a spherical tensor under 
rotation. § 

The commutators of the generators are most conveniently expressed in the spherical 
tensor coupled form (see appendix) 

[G(e)^^{/)](.) ^ _ ^_Y+f+9^gf^e f ,|^(,)^ (13) 

where we adopt the shorthand j = (2j + 1)^/^. The coefficient on the right hand side 
of (13) vanishes unless e + / + is odd. Consequently, the generators G^^ with g odd 
(^f = 1, 3, . . ., 2j — 1) close under commutation, forming the basis for the subalgebra 
SO(n). Finally, the generators G^\ which span the SO (3) algebra, are proportional to 
the physical angular momentum generators L^^^ — + l)(2j -|- l)]^^'^{a^ x a)^^^ for a 
single bosonic j-shell. 

For a single ferrnionic level, with degeneracy n = 2j + 1 (j half-integer), we have 
instead the chain U(n) D Sp(n) D ■ ■ ■ D SU(2) [see (2)]. The generators G^-f of U(n) 
again obey the commutation relations (13), and the generators with g odd {g = 1, 3, 
. . ., 2j) now span the Sp(n) algebra. The G^^^ are proportional to the physical angular 
momentum operators, now given by L^^^ = — [|i(j + l)(2j + l)]^/^(a^ x a)^x^ , closing as 
an SU(2) algebra. 

Proceeding now to the algebras involving both levels of the two-level system, let 
us reduce the complexity of the subscripts, relative to the generic multi-level notation 
^Im' by denoting the creation operators for the two levels by a]^(= a\ ,^) and fc^(= 
respectively, with angular momenta ja = ji and jb = j2, where a — —ja, —ja + 1; ■ ■ ■■, 

§ The convention ii"^ = (-)"+« aL"^ also arises in the literature. The relative sign between these 
conventions implies straightforward modifications A^"^ — >■ (— )2a to signs throughout the following 
results. 
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Table 4. Commutation relations for the generators of the U(ni + 712) algebra of the 

two-l('\"(4 hosonic s\"st(^ui (0 = +) or f( ritjionic s\-slom (0 = — ). in coupled form. 
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+ja and P — —jb, —jb + 1, . . ., +jb- The level degeneracies appearing in the algebra 
labels are rii — 2ja + 1 and ^2 = 2jb + 1. The algebra U(ni + ^2) is spanned by 

Gil^ = (at X a)(^) (^ = 0,l,...,2j„) 

Gi? = (at X (^=|i„-j6|,...,i„ + j6) .^4. 

G(f = (6t X a)(.'') (^=|i„-i,|,...,j„ + j5) ^ ^ 

Gif = (6t X (^ = 0,l,...,2j,). 

The commutation relations for these generators are listed in Table 4. They may 
all be obtained from the general bilinear commutation relation (A. 9). Notice the 
nearly identical commutation relations for the bosonic and fermionic realizations of 
the U(ni + 77,2) algebra, with sign differences indicated by the presence of the symbol 9 
in Table 4 (recall that 6 = + for the bosonic case and 6 = — for the fermionic case). 
Commutators not listed in Table 4, e.g., G\/J]'"^\ can be obtained from those given, 
by the coupled commutator symmetry relation (A. 4). The subalgebra Ui(ni) ® U2(n2) 
is obtained by simply omitting the "mixed" generators G^^^ and G^^J . 

The SO(ni + ^2) subalgebra, for the bosonic case, or Sp(ni + ^2) subalgebra, for 
the fermionic case, is then obtained by restricting (14) to the following generators: G^^a 
with g odd (^f = 1, 3, . . ., 2j„ - 1 or 2j„), G^f with g odd (^f = 1, 3, . . ., 2jb - 1 or 
2jh), and certain linear combinations of the form F*^^^ = rjgG^^j^ +^3^1^) coefficients 
of which are determined by the requirement of closure. Specifically, using the results 
of Table 4, it is found that closure is obtained if ig/rjg = o'o(— )^ for all g. That is, the 
relative sign between the terms must alternate, between generators F^^'^ with even and 
odd tensor rank g, but an overall sign parameter ctq may be chosen as either ±1. Thus, 
we obtain 

F^a) = ^j(at X 6)(f) + ao{-y{b^ x a)^^)]. (15) 
Note therefore that there are actually two distinct SO(ni + 712) subalgebras which may 
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Table 5. Commutation relations for the generators of the SO(ni + n2) algebra of the 
two-level bosonic system {9 = +) or the Sp(ni + n2) algebra of the fermionic system 
{6 = —), in coupled form. 



[i;(e)^i^(/)] 



•-Too 



^bb 
pif 

'^bb 



-2(-> 



^bb 





6l(_)ja+j5+/g/| 



ff <i f 9 \p(g) 

I ii, ja ja J 



Ja 3b 3h 



be included in (5), or two Sp(ni + 712) subalgebras in chain (6), distinguished by the 
relative sign ctq in the generators F'^^\ 

An overall arbitrary phase rig remains in the definition of F'^^\ If we choose F^^^ to 
be a "self-adjoint" tensor, in the sense that 

p{9)]^p{g)^ (16) 

i.e., Flf'^^ = {—)^^'^F^^, the commutation relations among the generators take on a 
simple form and, moreover, involve only real coefficients. Let ctq = (— with s = or 
1. Then a self-adjoint tensor F^^^ is obtained for the choice 

{1 a + b + s even , , 

i a + 6 + sodd. ^ ^ 

The commutation relations for the SO(ni + 712) or Sp(ni + 712) generators are listed in 
Table 5. 

The phase choice (17) for the two-level generator F^^^ also offers consistency with 
the s-b boson models, where F^-'*') plays an important role as a physical transition 
operator. For instance, in the IBM (jg — and jd — 2), the choice ao — + {i.e., 
s — 0) gives S0(6) generator F^^^ = (s^ x d)^^^ + (d^ x s)^^), which is the leading-order 
electric quadrupole operator [38]. The choice ctq = — {i.e., s = 1) instead yields the 
generator F*^^) = i[(s^ x d)^"^^ — {d'^ x s)*-^-*] of a distinct SO (6) subalgebra, denoted 
by SO (6) [53], which has been shown to be relevant to the decomposition of nuclear 
excitations into intrinsic and collective parts [54]. 

Finally, construction of the remaining subalgebras in (5) and (6) follows by 
application of the same principles. The SOi(ni) ® S02(n2) or Sp]^(ni) ® Sp2(n2) 
algebra is obtained by restriction to Cia"* and G^f^ with g odd, and S0i(3) S02(3) 
or SUi(2) (g) SU2(2) is obtained by further restriction to = 1. The combined angular 
momentum algebra, S0i2(3) or SUi2(2), then has generators L^^\ where 

L^'^ - elbaija + + l)]^/^G« + ^[b.O. + l)(2j, + 1)V^'G^^ ■ (18) 
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Table 6. Casimir operator eigenvalue formulas relevant to the two-level pairing model. Expressions are 
given first in terms of the generic Cartan irrep labels [A1A2 • • • Afe] and then specialized to the specific 
irreps arising for the two-level pairing system. 

Algebra (C2)[AiA2...] Irrep (C2) System 

SO(n) n = 2k + l 2X^{X^ + 2k + 1 - 2i) [v] 2v{v + n — 2) Bosonic single- level 

SO(n) n = 2k X^^Li 2A,(A, + 2A; - 2i) [v] 2v{v + n-2) Bosonic two-level 

Sp(n) n = 2k XlLi 2Aj(Aj + 2/c + 2 - 2i) {v} 2v{-v + n + 2) Fermionic 

U(n) Y^'l^^ + n + 1 - 2i) [N] N{N + n-l) Bosonic 

{N} N{-N + n + 1) Fermionic 



2v{ev + n-2e). (19) 



3.4- Casimir operators 

To exploit the symmetry properties of the two-level pairing model with respect to the 
subalgebras of U{ni +n2), it will be necessary (Sec. 5.1) to express the Hamiltonian in 
terms of the quadratic Casimir operators of the algebras in (5) and (6). Identification 
of the Casimir operator proceeds in two stages. First, a quadratic operator which 
commutes with the generators must be identified. This only defines the Casimir 
operator to within a normalization (and phase) factor. It is then desirable to choose the 
normalization such that the eigenvalues of the Casimir operator match the conventional 
eigenvalue formulas [43,55], given in terms of the Cartan highest-weight labels for the 
irrep in Table 6. For the symmetric irreps of SO(n) or antisymmetric irreps of Sp(n) 
arising in the two-level pairing problem, the eigenvalues can be expressed in terms of 
the single unified formula 

C2[S0(n)] 
C2[Sp(n)] 

Thus, as the second stage of defining the Casimir operator, the normalization is 
evaluated by explicitly considering the action of the operator on the one-body states, 
for which the irrep labels are known. 

For the single-level algebra, SO(n) or Sp(n), the operator G oG, defined by|| 

GoG^-Y^ 9[G^'^ X G^V, (20) 

g odd 

commutes with all the generators, i.e., G^^^ with g odd [24,31,46]. The result follows 
from the general theory of Casimir operators for an algebra, and it may be verified 
by explicitly evaluating the commutator [G'^^\G o G]^^\ using commutator (13) and 
product rule (A. 6). In terms of the conventional spherical tensor scalar product, defined 
by A(^) • S(^) = (-)^c(A(^) X S(^))J°\ this operator is G o G = EsoddG"^^^ ' G^'^^- The 
eigenvalue olGoG acting on the one-body state \a) may easily be evaluated by Wick's 

II The generators of SO(n) or Sp(n) together transform as an SO(n) or Sp(n) tensor. Therefore, 

the circle in the notation G o G is meant to represent a scalar product with respect to SO(n) or 
Sp(n), following Ref. [56]. In (22) the notation is generalized to represent an SOi(ni) ® 802(712) or 
Spi(ni) ® Sp2(n2) scalar. 
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4GoG, (21) 



29FoF + AGaa O Gaa + ^G^b o G^b (23) 



theorem in coupled form [57], using the commutator results described in the appendix. 
Comparison with the eigenvalue formula (19) gives normalization 

C2[S0(n)] 
C2[Sp(n)] 

covering both the bosonic and fermionic cases. 

Proceeding to the two-level problem, for the Casimir operator of SO(ni + or 
Sp(ni + n2), we start from the Casimir operators AGaa ° Gaa and 46*55 ° ^66 of each 
single-level subalgebra, as defined in (20). Although each of these operators commutes 
with each of the generators of SOi(ni) Cg) S02(n2) or Sp]^(ni) (X> Sp2{''^2), they do not 
commute with the two- level generators F^^\ We therefore introduce the operator 

FoF = ^^[F(^)xF(^)](°\ (22) 

9 

or, equivalently, FoF = '^gi—^F^^'' ■ F^^\ with g = \ja — jb\, ■ ■ ■ ,ja + jb- This quantity 
is invariant with respect to SOi(ni) ® 802(^2) or Sp;^(ni) (X> Sp2(n2). Moreover, the 
combination 

C2[SO(ni+n2)] 
C2[Sp(ni + n2)] 

commutes with all the generators of SO(ni -\-n2) or Sp(ni +n2), as seen by application 
of the commutators in Table 5 and the product rule (A. 6). That this combination 
of operators also has the correct normalization to match the eigenvalue formula for 
C2[S0(n)] {n even) or C2[Sp(n)] (Table 6) may be verified by explicitly calculating the 
one-body expectation value {a\G2\a) or {b\C2\b). 

Returning to the example of the IBM SO (6) D SO (5) chain, the Casimir 
operator (21) becomes 

C2[SO(5)] = Aid'' X J)« • {d^ X dy^) + A{d^ X df^ ■ {d^ X df\ (24) 

for the single level consisting of the quadrupole boson d^'^\ Then, for the S0(6) algebra 
of the two-level s-d system, 

C2[SO(6)] =2[(stxJ)(2)+(dtxs)(2)].[(stxJ)(2)+(rftx5)(2)]+C2[SO(5)],(25) 

consistent with the usual result [22]. 

Similar results apply to the quadratic Casimir operators of the unitary algebras 
in (5) and (6). For the single-level algebra, the linear invariant of U(n) is simply the 
occupation number operator N — clI^cl^, or N — 9jG^q \ The quadratic invariant is 
given by 

C2[^{n)] = Y^g{-y[G^9) X G'(^)](o), (26) 
9 

or C2[U(n)] = YligG''^'' ■ G^^\ However, for the bosonic realization of U(n), only 
symmetric irreps arise, and, for the fermionic realization of U(n), only antisymmetric 
irreps arise, with eigenvalues as given in Table 6. Therefore, in either situation, it is 
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found that the quadratic invariant is simply a function of the linear invariant and can 
be expressed as 

C2[U{n)]^N{eN + n-e). (27) 

Likewise, the Casimir operator of the two-level system's algebra U(ni + ^2) may be 
expressed as 

C2[U(ni + 722)] = 2eNiN2 + n2Ni + miVa + C2[Ui(ni)] + C2[U2(n2)]. (28) 

This result is obtained by comparison of the eigenvalues for C2[U(ni + 722)] with those 
for C2[U(ni)] and C2[U(n2)], together with the additivity of the number operators 
{N^N^ + N2). 



4. Quasispin algebra 

First, we note that a set of three operators S+, S^, and obeying the commutation 
relations 

[^0,^+]=+^+ [^0,^-] = -^- [S+,S^]^-2eSo, (29) 

and obeying the unitarity conditions 5"^ = and Sq — Sq, span a unitary realization 
either of the algebra SU(1,1), for ^ = +, or of the algebra SU(2), for 9 ^ -. The 
SU(1, 1) or SU(2) invariant operator is given by 

Sq — ^9{S+S- + S-S+) , , 

^So(So-i)-es+S-. ^"'''^ 

For an irrep of SU(1, 1), this operator takes on eigenvalues S{S — 1), and the possible 

eigenvalues of 5*0 are given byM = S, -5-1-1, S'-|-2, For the "true" group 

representations of SU(1,1), S must be integer or half-integer, but description of the 
bosonic system in the quasi-spin formalism as considered below requires the projective 
representations with S = |, |, |, for reasons described in Ref. [34]. As usual, for 
SU(2), takes on eigenvalues S{S + 1), with 5 = 0, |, 1, . . ., and the eigenvalues of 
5*0 arc given by M = —S, . . ., +S — 1, -l-^. 

Now, consider a system consisting of one or more j-shells of angular momentum jk 
{k — 1, 2, . . .). Regardless of whether the operators a^^ for the kth level are bosonic or 
fermionic, a quasispin algebra is defined following the prescription of Sec. 2. The scalar 
pair creation operator, scalar pair annihilation operator, and an operator simply related 
to the number operator for this level form a closed set under commutation. Specifically, 
let 



SkQ — \ ^nX'^\nfikm + ^OfemO^ 



km) 



These define cither an SU(1, 1) quasispin algebra for bosons [34], which we denote by 
5Wfc(l, 1), or an SU(2) quasispin algebra for fermions [25,26], which we denote by 
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SUki'^). I" terms of spherical tensor coupled products, these operators (31) may be 
represented as 

Sk+ = iife(4 X «I)o°^ 

Sk- - \9U^k X ~akt^ (32) 
Sko^\3ke[{alxakt'^ + {akxalt\ 

The operator Sko is related to the number operator Nk — ^Im^fe"i' which may be 
expressed in spherical tensor form as Nk — ife(a[ x afe)Q°\ by 

Sko^UNk + eQk), (33) 

where the constant flk = \{'^jk + 1) is the pair degeneracy of the level k. The quasispin 
Ski moreover, is related to the seniority quantum number Vk by the duality relation for 
a single j-shell. 

It is therefore possible to interconvert between quasispin quantum numbers Sk and 
Mki for each level, and occupation-seniority quantum numbers Nk and Vk, according to 

Sk^\{^k + Ovk) , . 

Mk^l{Nk + enk). ^ ^ 

Since the lowest weight state for a given quasispin (i.e., with Mk — 9Sk) is destroyed 
by the pair annihilation operator, and since this state contains Nk — Vk particles, Vk 
may be interpreted as the number of unpaired particles, either bosons or fermions. The 
seniority Vk takes on values Vk = 0, 1, 2, . . ., subject to the constraint Vk < Nk for 
bosons or Vk < min(A^fc, 2Qk ~ Nk) for fermions, by the M contents of the irreducible 
representations noted above. ^ Note, therefore, in the fermionic case, that Vk < flk, with 
the maximum value occuring for half filling (A^^ = flk)- 

A quasispin algebra for the two-level system — which we denote by SUi2{l, 1) for 
the bosonic case or 5Wi2(2) for the fermionic case — is spanned by the sum generators 

5+ = Si+ + 0-5*2+ = Si^ + aS2- So = Sio + S2o- (35) 

A quasispin algebra is obtained with either choice of sign cr = ± in the ladder operators. 
This algebra defines a total quasispin quantum number S which is dual to the two-level 
algebra seniority quantum number v, by a relation of the same form as (34), namely, 

s = l{Q + ev) 

M='^iN + dn), ^^^^ 

where N and Q are defined above as the sums of the single-level values. The allowed 
values for the total quasispin S are given for 5Wi2(l, 1) by 5* > Si + S2 {i.e., 

^ As noted in Sec. 3.1, the case of a j = bosonic level is anomalous. There is no orthogonal algebra 
dual to the quasispin algebra, and thus no seniority quantum number, but the label Vk may still be 
defined from the quasispin via (34). The squared quasispin operator is identically = — ^ for such a 
level, as described in Refs. [34,58]. Thus Sk = j 01 |, and hence Vk = ov 1. Since M — S must be 
integral, it follows that Wfc = for Nk even and Vk = 1 for Nk odd, i.e., Vk = Nk mod 2. The natural 
interpretation of this value is that particles in a j = level axe automatically paired to zero angular 
momentum, except for the one unpartnered particle when the occupation is odd. 
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5*1 ® 5*2 — )■ + S2, + S'2 + 1, . . .) and for SUi2{2) by the familiar triangle inequality 
\Si -S2\<S<Si + S2 {i.e., Si® \Si - ^2!, S1 + S2-I, Si + S2). 

When reexpressed in terms of the seniority labels Vi, f 2, and v, the SU{1, 1) coupling 
rule is equivalent to the SO(ni + ^2) — )■ SOi(ni) ® 802(^2) branching rule (8), and the 
SU{2) coupling rule is equivalent to the Sp(ni + 712) Spi(ni) ® Sp2(n2) branching 
rule (10). Similarly, the SU{1,1) M content rule is equivalent to the U(n) — )■ SO(n) 
branching rule (7), and the SU{2) M content rule is equivalent to the U(n) — )■ Sp(n) 
branching rule (9). 



5. Hamiltonian relations 

5.1. Dynamical symmetries 

Before considering the pairing Hamiltonian in particular, let us consider the Hamiltonian 
defined by the Casimir operators of the unitary subalgebra chains. A dynamical 
symmetry [43] arises when the Hamiltonian is constructed in terms of the Casimir 
operators of a single chain of subalgebras. The eigenstates thus reduce the subalgebra 
chain, i.e., constitute irreps of the subalgebras. More generally, especially when 
considering phase transitions, it is useful to construct the Hamiltonian from terms 
consisting of Casimir operators from multiple, parallel chains, here the upper and lower 
chains of (5) or (6), as 

C2[SO(ni + n2)] 
C2[Sp(ni + n2)] 



H = aN + biNi + b2N2 + b 



where higher-order invariants may also be included 

The upper chain in (5) or (6) defines an SO(ni + 71,2) or Sp(ni + 77,2) dynamical 
symmetry, and the lower chain defines a Ui(77.i) U2(7i2) dynamical symmetry. The 
Ui(7ii) 112(772) dynamical symmetry is obtained for the Hamiltonian (37) with 6 = 0, 
i.e., 

H = aN + biNi + b2N2 

( C2[S0,{n,)] 1 r C2[S02(772)] 1 , J2 , , j2 , .2 .30^ 

C2[Sp,(770] C2[Sp2{n2)] ! + ^^-^^ + ^^-^^ + ■ ^^^^ 

The eigenstates | A^A^iA^2^i^2 ■ ■ ■ J1J2J) have definite occupation numbers for each of the 
levels and have energy eigenvalues 

E = aN + biNi + b2N2 + 2ciVi{evi + 77i - 29) + 2c27;2(^t'2 + r72 - 2d) 

+diJi{Ji + 1) + ^2^2(^2 + 1) + eJ{J + 1). (39) 

The S0(7ii + n2) or Sp(ni + 772) dynamical symmetry is obtained for the 
Hamiltonian (37) with 61 = 62 = 0, i.e., 

C2[SO(77i+ 772)] 



H^aN + b 



C2[Sp(77i +772)] 
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Figure 1. Energy diagrams for the bosonic SO(ni -|- 712) and fermionic Sp(ni -|- 
712) dynamical symmetries of the two- level pairing model. The degeneracy within 
irreducible representations of the two-level algebra (brackets labeled by v) is split 
according to the single-level algebra irreducible representation labels {viV2). (a) Energy 
levels for the bosonic SO (6) D SO (3) (g) SO (3) dynamical symmetry Hamiltonian 
(ii = J2 = !)• (b) Energy levels for the fermionic Sp(8) D Sp(4) (g) Sp(4) dynamical 
symmetry Hamiltonian (ji = j2 = §)• The Hamiltonian in each case is chosen as 
H = —OF o F, i.e., b = and Ci = C2 = +^ in (40). The total occupation in both 
cases is N = A, which for the fermionic example gives half filling. 



with eigenstates |A^f f 1^2 ■ ■ ■ J1J2J) and energy eigenvalues 

E^aN + 2bv{ev + ni + ^2 - 2d) + 2ciVi(dvi + ni - 2d) + 2c2V2{9v2 + n2 - 26) 

+diJi{Ji + 1) + 4^2(^2 + 1) + eJ( J + 1). (41) 

The energy spectrum for the dynamical symmetry follows from the branching rules of 
Sec. 3.2. Example level energy diagrams are shown for a bosonic system (ni = n2 = 3) 
in Fig. 1(a) and for a fermionic system of similar degeneracies {rii — 712 — 4:) in Fig. 1(b). 



5.2. Pairing Hamiltonian 

The pairing Hamiltonian for a generic multi-level system, consisting of levels of angular 
momentum ji, J2, • • •) is given by 

^ = I]£ife«L«ifcm+i G'ifc'fe(-)"='""'4'm'4',-m'(-)''""«fc,-m«ikm,(42) 

km k'k 
m'm 

where the summation indices k and k' run over the single-particle levels (/c = l, 2, ...), 
and m and m' run over their substates (m = —jk, —jk + 1, ■ ■ +jk)- The creation 
operators aj[.^ and annihilation operators Ukm obey either canonical commutation 
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relations (bosons) or anticommutation relations (fermions). The first term represents 
the one-body energy contribution for each level, and the products in the second term are 
creation or annihilation operators for pairs involving time-reversed partner substates. 
Note that Gk'k = G\j^, by Hermiticity of H , and that we take all coefficients to be real. 
That this Hamiltononian is specifically constructed from angular momentum zero pair 
operators is seen by rewriting it in terms of angular-momentum coupled products, as 

H ^^EkNk + \e^Gk'kJk'Ual' X 4')o\^k X au)f. (43) 

k k'k 

The Hamiltonian (42) is integrable and may be solved using a generalized 
Gaudin algebra [59] or, under certain conditions, Bethe ansatz [60, 61] methods. Two 
limiting cases deserve special mention, since they are characterized by dynamical 
symmetries (Sec. 5.1), as described below, and are solvable by more elementary methods. 
(1) Trivially, when all Gk'k = 0, the problem reduces to that of a system of noninteracting 
particles (weak-coupling limit), with eigenstates characterized by occupations numbers 
Nk- (2) When all Ek = 0, physically corresponding to the situation in which the level 
energy difference is negligible relative to the pairing strength (strong-coupling limit), 
and if all Gk'k are equal to within possible phase factors (uniform pairing), the problem 
is immediately solvable by the use of quasispin (Sec. 5.3). 



5.3. Quasispin Hamiltonian 

The generic multi- level pairing Hamiltonian (42) can be expressed entirely in terms of 
the quasispin generators, by comparison with (32), as 

H = ^ Ski'^Sko - 9Vtk) + ^ Gk'kSk'+Sk-- (44) 

k k'k 

This Hamiltonian therefore conserves the quasispin Sk or, equivalently, the seniority Vk, 
associated with each level. The Hamiltonian is also number conserving, so the total z- 
projection 5*0 = J2k ^ko — ^{N + 9Q) is conserved, where N — J2k and Q = J2k ^k- 
However, the individual Sko are not in general conserved, unless the levels completely 
decouple, with Gk'k = for all k' ^ k. 

Note that numerical diagonalization is straightforward in the weak-coupling basis, 
consisting of states |A''if ■ ■ ■) of good seniority and occupation for each level. The 
action of the Hamiltonian (44) on these states follows from the known action of the 
quasispin ladder operators, S±\SM) = [-e{S ± 9M){S ^OM - e)Yl^\S{M ± 1)), once 
the quantum numbers are translated via (34). Specifically, 

Sk^\ ■■■NkVk---) = \[e{Nk -Vk + 2){Nk + Vk + 2enk)Y'^\ ■■■{Nk + 2)vk ■ ■ ■) , . 
Sk-\--- NkVk ■■■) = l[e{Nk - Vk){Nk + Vk + 2dnk - 2)]V2| ■ ■ ■ (iV^ - 2)vk ■■■). ^ ^ 

Consequently, the matrix elements for the diagonal pairing terms are 

(• • • NkVk ■ ■ ■ \Sk+Sk-\ ■ ■ ■ NkVk ■■■) = \e[Nk{Nk + 2enk - 2) - Vk{vk + 2enk - 2)] (46) 
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and for the off-diagonal terms are 

(• • • (Nk' + 2)vk' ■■■{Nk- 2)vk ■ ■ ■ \Sk'+Sk-\ ■ ■ ■ Nk'Vk' ■ ■ ■ N^Vk ■ ■ ■) 
= \[{Nk' - Vk' + 2)(Nk' + Vk' + 2enk') 

x{Nk - Vk)iNk + Vk + 2enk - 2)]V2, (47) 

as noted for the fermionic case in, e.g., Refs. [10,40]. 

Returning to the two-level problem, the pairing Hamiltonian in quasispin notation 

is 

l+Si- + £2(25'20 — ^^^2) + G22S2+S2- 

+Gi2{Si+S2- + S2+S,-). (48) 

The two-level pairing Hamiltonian has two dynamical symmetries [62] defined with 
respect to the quasispin algebras, corresponding to either the upper or lower subalgebra 
chains in 

<SWl2(l,l) 

SUi(l, 1) (g) 5^2(1, " } dUi2(1) (49) 



Vl V2 



Ui(l)®U2(l) I N 



for the bosonic case or 

SUi{2) ® SU2{2)D { TTn^ITTn^ } ^^i2{^) (50) 




for the fermionic case, with conserved quantum numbers as indicated. Here the algebra 
Ufc(l) is the trivial Abelian algebra spanned by Sko, and Ui2(l) is spanned by the their 
sum 5*0, as defined in (35). The occupation-seniority labels {v and N) are indicated, 
rather than the quasispin labls {S and M), for a closer connection to the physical 
problem and easier comparison with the dual algebra's dynamical symmetries. 

The dynamical symmetry Hamiltonian for the upper subalgebra chain is the strong 
coupling limit (i.e., £1 = £2 = 0) of the two- level pairing Hamiltonian, with uniform 
pairing strength, as defined in Sec. 5.2. Specifically, let Gn — aGu — (7G21 — G22 = C, 
for either sign o" = ±. Then the Hamiltonian is given by 

H^GS+S^, (51) 

where S± are the sum-quasispin ladder operators of (35), defined in terms of the same 
sign a. Since S+S- — 9[So{So — 1) — S^j, by (30), the strong-coupling Hamiltonian 
conserves the total quasispin S (or seniority v), as well as the projection quantum 
number M (or occupation N), and has eigenvalues 

= \e[N{N + 2en - 2) - v{v + 2en - 2)], (52) 

as expressed in terms of the occupation-seniority labels. The eigenstates are identical to 
those of the SO(ni -|- 712) or Sp(ni -|- ^2) dynamical symmetry Hamiltonian (40)."*" The 

+ More precisely, the quasispin Hamiltonian (51) has a higher degeneracy than the SO(ni -|- 712) or 
Sp(ni +^2) Hamiltonian (40), but the eigenstates can be chosen from within each degenerate subspace 
to match those of (40), i.e., of good Ji, J2, and J. 
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specific relationship between the Hamiltonians is determined below in Sec. 5.4. 

The dynamical symmetry Hamiltonian for the lower subalgebra chain is the weak- 
coupling limit of the pairing Hamiltonian {Gk'k = 0), as defined in Sec. 5.2. The 
dynamical symmetry eigenstates are simply the level occupation eigenstates of good Ni 
and A''2, as for the Ui(ni) ® 112(^2) dynamical symmetry Hamiltonian (38), i.e., the 
weak-coupling basis states considered above. 

The full two-level pairing Hamiltonian (48) can be expressed entirely in terms of 
the invariant operators of algebras appearing in the upper and lower chains. Specifically, 

H = £i(25io - en^) + (Gn - ^G'i2)^[^io(^io - 1) - S^] 

+£2(2^20 - ^^^2) + (^22 - CTGu)e[S2o{S20 - 1) - S^] 

+Gi2ae[So{So - I) - S% (53) 



5.4. Duality relations for the Hamiltonian 

The eigenstates for the dynamical symmetries of the two algebraic frameworks — 
number-conserving unitary and numbcr-nonconserving quasispin — are identical, that 
is, the irreps which reduce the unitary algebra chains (5) and (6) reduce the quasispin 
algebra chains (49) and (50) as well, and the labels for the chains are connected through 
the duahty relations. The pairing Hamiltonian is defined in Sec. 5.2 [see (42)] in 
terms of certain combinations of operators which represent scalar pair creation, scalar 
pair annihilation, and number operators. These are noted in Sec. 5.3 [see (44)] to be 
essentially the quasispin generators, and the Hamiltonian can also therefore be expressed 
directly in terms of the quasispin invariants [see (53)]. However, the pairing Hamiltonian 
can just as well be expressed in terms of the Casimir operators of subalgcbras of 
U(ni-|-n2), as a special case of the Hamiltonian of Sec. 5.1 [see (37)]. That such a relation 
exists is implied by the duality of irreps, but it is explicitly obtained by appropriate 
recoupling and reordering of the bosonic or fermionic creation and annihilation operators 
in this section. 

For a single j-shell, recoupling and commutation of creation operators* yields 
(at X at)(°)(a x a)^^^ 

^ _(at X a)i°) + Y,g[{a^ X a)(^) x (at x d)^^^^i'^ 



(54) 

Thus, the relation between quasispin and Casimir Hamiltonians for a single level is 

,S,S.^-0N^C^i^-l[^^], (55) 

by comparison with the explicit realizations of the various operators, namely, N from 
Sec. 3.4, C2[S0(n)] or C2[Sp(n)] from (21), C2[U(n)] from (26), and S± from (32). 

* The product of a pair creation operator and a pair annihilation operator is related to the product of 
spherical-tensor one-body operators by, e.g., identity (25a) of Ref. [57]. 
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For the two-level system, which has quasispin generators given by (35), the product 
S^S^ involves both single-level terms (5'i+5'i_ and S2+S2-) and "cross terms" (crS'i+S'2_ 
and aS2+Si^), which destroy a pair in one level and create a pair in the other. 
Recoupling and commutation of the mixed terms yields 

(at X at)(°)(6 X b)^^^ + {b^ x b^)^^\a x 5)^°^ 

= ^ \ . ' [FoF- 2NaN, - 9{f,Na + fM] . (56) 

Jajb 

The one-level terms and mixed terms of S+S^ may thus be combined to give an 
expression involving the SO(ni -|- 712) or Sp(ni -|- 712) Casimir operator and U(ni -|- 712) 
invariants, if and only if the sign a arising in the definition of the sum quasispin algebra 
and the sign ao entering into the definition of F^^^ are related by 

^ = -g(-ya+jt_ (57) 
(J 

We again have an expression of the same form as (55), 



That the expression is of this form is to be expected from the general nature of the 
duality, which indeed makes no assumption (Sec. 2) as to whether the single-particle 
states are considered to be arranged into a single j-shell, as for (55), or two j-shells, 
as here. Prom a practical standpoint, what is most useful is that the relation can 
now be expressed explicitly in terms of the two- level system operators given by (23), 
(28), and (35) as spherical-tensor products of creation and annihilation operators, with 
well-defined phases ao and a. 

The two-level operator correspondence (58) relates the 5^12(1, 1) or SUui^) 
dynamical symmetry Hamiltonian in the quasispin scheme (51), i.e., strong coupling 
with uniform pairing, to the SO(7ii-|-7i2) or Sp(7ii-|-7i2) dynamical symmetry Hamiltonian 
in the unitary algebra scheme. ji Furthermore, taken in conjunction with the single-level 
operator correspondence (55), it allows the full two-level pairing Hamiltonian (not just at 
the dynamical symmetry limit) to be expressed in terms of Casimir operators of algebras 
appearing in the two parallel subalgebra chains of U{ni -\- 712)- Starting from (48), one 
obtains 

//=(£!- leGn)N, + |(Gn - aG,2)C2[\JM] - l(Gn - c7Gu)\ ^/t^f'^^, 
+(£2 - 10G22)N2 + 1(^22 - <7G^2)C2[U2{n2)] - 1(^22 - aGi2)| ^^[Sp^rij] 

H^O.^C2mn.^n2)^-l^G.^^^^^^^^^ (59) 

jj The analog of (58) for the IBM was exploited in Ref. [38] to establish the properties of the IBM 
S0(6) dynamical symmetry eigenstates. With = the sign condition (57) gives ctq = —(^ (recall 
ja = 0, and jb = 2). Thus, the pairing operator S+S- for the quasispin defined with negative sign 
(cr = — ) relates to the Casimir operator of the "physical" S0(6) algebra ((Tq = +), and the quasispin 
algebra defined with positive relative sign is instead dual to the S0(6) algebra (Sec. 3.3). 
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For this relation to be valid, the phases Uo and a, used in defining the generators for 
orthogonal or symplectic algebra and two-level quasispin algebra, respectively, must be 
related by (57). 



5.5. Multipole Hamiltonian 

The Hamiltonian for spectroscopic studies of the s-h boson models is commonly 
expressed in terms of a "multipole" term of the form [{s'^ x 6)(^) + (6''' x s)^-^)] ■ \[s^ x 
6)(-^) + (6t X s)(^)], where L = jb [63]. For instance, the customary IBM U(5)-SO(6) 
quadrupole Hamiltonian [64, 65] is 

The U(5) limit is obtained for ^ = and the S0(6) limit for { = 1. The operator FoF, 
appearing as the "cross term" in C2 [SO (711+712)] or C2[Sp(7ii+7i2)] [see (23)], generalizes 
the multipole term to the generic two-level model. In contrast, the Hamiltonian for 
spectroscopic studies of the fermionic system is commonly expressed in pairing or 
quasispin form [5, 20]. Thus, we seek to relate these distinct — pairing and multipole — 
forms of the Hamiltonian. 

Recall that the operators considered thus far in connection with the strong-coupling 
limit — S^, S+S-, and C2[SO(72i +712)] or C2[Sp(7ii +n2)] — differ only in normalization 
(or sign) and by addition of a function of A^, the conserved total occupation number. 
Therefore, the eigenstates are identical and the eigenvalues differ only by a rescaling and 
a constant offset. However, the operator FoF appearing in the multipole Hamiltonian 
differs from these [again, see (23)] by terms proportional to C2[SOi(7ii)] and C2[S02(7i2)], 
in the bosonic case, or C2[Spi(7ii)] and C2[Sp2(7i2)], in the fermionic case. The eigenstates 
arc therefore again the same as for S^, S+S^, and C2[S0{ni + 712)] or C2[Sp(7ii + n2)], 
but eigenvalues are no longer degenerate for states sharing the same value of v. They 
are rather now split by Vi and 7;2, as illustrated in Fig. 1. 

For the explicit relationships among these operators, observe that, by (58), 

and, in terms of (23), 

C2[Sp(7li+7l2)] J I C2[Sp,{m)] C2[SP2(712)] j'^ ' 

Therefore, the pairing (S+S^) and multipole {F o F) forms of the Hamiltonian, those 
most frequently encountered in applications, are related by 



-0FoF^ AS^S^ 



C2[U(7li+7l2)] -^A^ 



ir C2[SOi(770] 1 If C2[S02(772)] . . . 

^21 C2[Spi(7Zi)] /^2l C2[Sp2(n2)] ^' ^ ^ 
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where 

simply contributes a c-number shift to the eigenvalue spectrum, without affecting the 
eigenfunctions. 

A positive coefficient (G > 0) for S+S- gives a positive pair energy, i.e., repulsive 
pairing, in both bosonic and fermionic cases, as may be seen from (52) with N = 1 and 
V = 0. The sign of the pairing interaction is of special interest in comparing the bosonic 
and fermionic two-level pairing models, since it should be noted (Sec. 6) that the system 
undergoes a quantum phase transition for repulsive {G > 0) pairing interaction in the 
bosonic case and attractive pairing interaction {G < 0) in the fermionic case. Thus, it 
is essential to note that repulsive pairing is obtained for a negative coefficient on F o F 
in the bosonic case and a positive coefficient on F o F in the fermionic case, i.e., for 
H = —OF o F, or vice versa for attractive pairing. Therefore, for repulsive pairing, a 
Hamiltonian 

HpF = ^^^2 - ej^^F o F (65) 

is the natural generalization of the multipole form for the transitional Hamiltonian (60) 
to generic two-level pairing models, as considered in Sec. 6. 

The last two in terms in (63), involving the Casimir operators of SOi(ni) and 
S02(n2) or Spi(ni) and Sp2(n2), contribute a common shift to the energy eigenvalues for 
each subspace of states characterized by a given pair of values of the conserved (^1^2) 
quantum numbers, without affecting the eigenfunctions, i.e., these terms serve only to 
displace the different (f 1^2) subspaces relative to each other. If only {V1V2) = (00) states 
arc considered, the SOi(ni) and 802(^^2) or Spx(ni) and Sp2(n2) terms have no effect at 
all. They will therefore not be considered further. 

Now to consider the spectra, the strong-coupling Hamiltonian operator AS+S^ — we 
include the factor of 4 arising in (63) for convenience — has eigenvalues given by (52), 
obtained with < v < N for bosonic pairing or with < v < min(A'", 2Q — N) for 
fermionic pairing, where only even values of v arise for N even, or odd values of v for 
N odd [see (7) and (9)]. Thus, taking N even, the eigenvalues span the range 

(AS+S-)^ ,...,N{N + 2n-2), (66) 

for bosonic pairing, or 

(AS+S-)^ ^ or A(N -n) ,..., N(2n- N + 2\ , (67) 

^ = ^ v = 2n-N v=0 

for fermionic pairing, in which case < < 2Q. (If N is odd, the sequences above would 
end instead with v = 1 rather than v = 0, but the large- iV dependence of the highest 
eigenvalue on A^^ is not changed.) The range of eigenvalues therefore depends upon the 
total occupation or "ffiling" A^ of the two- level system as sketched in Fig. 2(a) for the 
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Figure 2. Eigenvalues of the pairing interaction term 45"+ S*-, which determines the 
energy spectrum of the two-level pairing model in the strong coupling limit, shown for 
the (a) bosonic and (b) fermionic systems, as a function of filling N. The eigenvalues 
are given by (66) and (67), respectively. The axes arc labeled gcnerically, to indicate 
the asymptotic (large-fJ) dependences discussed in the text, but the specific points 
shown for illustration are calculated for f2 = 50. 

bosonic system and Fig. 2(b) for the fermionic system. The asymptotic dependences for 
large degeneracy (Q » 1) are indicated. Specifically, at a fiUing approximately equal to 
half the total degeneracy, i.e., N !=i note that the bosonic eigenvalues span a range 
~ 3Q^(Ri 3A^^), while the fermionic eigenvalues for the same fiUing and degeneracy only 
span a range of ~ N'^). If, instead, the limit of large occupation is taken at fixed 

degeneracy (Q ^ N) in the bosonic case, the range of eigenvalues ~ N"^ is the same as 
for a fermionic pairing model of the same N but at half filling (which is obtained for a 
correspondingly larger degeneracy Q = N). 

Note also that, for repulsive pairing, the bosonic ground state (for which v = N) 
has zero eigenvalue. In contrast, the fermionic ground state (for which v — below half 
filling and v — 2D, — N past half filling) has an eigenvalue which grows linearly with 
N past half filling. The nonzero ground state pairing energy for the fermionic system 
may be understood since, past half filling, Pauli exclusion enforces the existence of some 
particles in time-reversal conjugate orbits and hence some probability for pairs coupled 
to zero angular momentum. 

For the multipole form —6F o F of the pairing interaction operator, the spectrum 
is shifted downward by an iV-dependent offset relative to that of 45'+S'_ [see (63)]. The 
highest eigenvalue (obtained for i> = 0) is always zero. The asymptotic form of the 
ground state eigenvalue is, alternatively, ~ —N"^ for fermionic half fiUing {1 — N), 
~ —SN'^ for bosonic "half fiUing" (1 <^ Q fa A?"), and again ~ —N"^ for the bosonic 
system at larger boson number {Q <^ N). 
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6. Transitional Hamiltonian 

A second- order quantum state phase transition occurs between the weak-couphng 
and strong-couphng hmits for the two-level pairing models. Specificially, for the 
bosonic system it occurs with repulsive pairing interaction, and for the fermionic 
system it occurs with attractive interaction. The quantum phase transition is apparent 
numerically from calculations for finite N and from semiclassical treatments of the large- 
N limit. The present duality relations (Sec. 5) immediately help clarify comparison of 
numerical eigenvalue spectra across the transition but are also intended to facilitate the 
construction of coherent states for the semiclassical treatment. 

The simplest semiclassical "geometry" for the two-level pairing model is obtained 
from the quasispin algebraic structure, most simply by replacing the quasispin operators 
with classical angular momentum vectors, which maps the pairing model onto an 
essentially one-dimensional coordinate space. This approach has been applied in both 
the bosonic s-b models and fermionic two-level pairing model with equal degeneracies [5, 
8,42,66,67]. In both these circumstances, the quantum phase transition is found to 
occur, in the large- A'^ limit, at N\G\/e = 1. For the s-b models, a higher- dimensional 
and richer classical geometry (see Ref. [42]) has been established through the use of 
U(n2 + 1)/U(n2) coherent states [42,68-71]. An extension of this treatment to U(ni -|- 
^2)/[U(ni) ® U(n2)] for generic two-lcvcl pairing models might profitably be obtained 
using the exphcit construction of generators for the SO(ni -|- 712) D SO(ni) ® S0(n2) 
and Sp(ni -|- n2) D Sp(ni) ® Sp(n2) chains considered in Sec. 3. 

However, at present, we confine ourselves to laying the groundwork for more 
detailed further work, allowing for the most general choice of level degeneracies and 
more uniformly treating the bosonic and fermionic cases. A pairing Hamiltonian 

^pair = ^^iV2 + ^i|^+5_, (68) 

may be defined with opposite signs 9 of the pairing term for the bosonic and fermionic 
cases, so that the quantum phase transition is obtained in either case. This Hamiltonian 
yields the weak-coupling limit for ^ = 0, the strong-coupling limit at ^ = 1, and the 
critical interaction strength N\G\/s = 1 at ^ = 1/5. Scaling of the one-body term 
by A^"^ and of the two-body term by by N~'^ ensures that the critical point remains 
fixed at the finite value ^ = l/5asA^— >oo. However, for this Hamiltonian, a grossly 
different "envelope" to the eigenvalue spectrum (i.e., the range of eigenvalues, obtained 
as a function of the control parameter ^) is found in the bosonic and fermionic cases 
(see Fig. 3). 

To facilitate direct comparison of the bosonic and fermionic quantum phase 
transitions, it is helpful to instead construct a Hamiltonian for which the ground state 
energy follows the same trajectory as a function of ^ in the large A^ limit, and the 
eigenvalues span the same range at each of the hmits, namely, [0, 1] for ^ = and [—1, 0] 
for ^ = 1. By the results of Sec. 5, this is accomplished by choosing, for repulsive pairing 
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Figure 3. Lowest and highest eigenvalues, within the {V1V2) = (00) subspace, for the 
bosonic two-level pairing model with repulsive pairing interaction, in the usual form as 
given by -ffpair in (68) (solid curves) and after correction by the U(ni-|-n2) Casimir offset 
as given by in (69) (dashed curves). These are plotted as functions of the control 
parameter ^ between the weak-coupling and strong-coupling limits. The eigenvalues 
for the fermionic two-level pairing model with attractive pairing interaction, as given 
by -ffpair or equivalently by H- in (70), are shown for comparison (dotted curves). All 
calculations are for = 50, with m = 722 = 5 in the bosonic case and ni = 722 = 50 
in the fermionic case (see Fig. 4 for more detailed eigenvalue spectra). 



45 , - 



(69) 



interaction, the Hamihonian 

N{N + 2VL-2) 
N{2VL -N + 2) 

and, for attractive pairing interaction, the usual Hamiltonian 

d _ f) 4£ 

H--^-^N,-f^S^S^. (70) 

The c-number offset included in the definition of which arises as C2 [U(ni 4-77,2)] — 6* A^, 
is included to achieve the same range of eigenvalues in the strong coupling limit, as well 
as a similar evolution of ground state energy across the transition (Fig. 3), in the large- 
N limit, thereby facilitating comparison of the bosonic (repulsive pairing) and fermionic 
(attractive pairing) quantum phase transition. With inclusion of this offset, if+ is 
equivalent to the generalized multipole transitional Hamiltonian (65) when acting on 
the {V1V2) = (00) subspace of any two- level pairing model, as may be seen from (63). 
In particular, for the s-b models, inclusion of this offset makes identical to the 
conventional multipole form (60) of the transitional Hamiltonian, when acting on the 
{V1V2) = (00) subspace. 

The evolution of the eigenvalue spectrum across the transition between weak 
coupling and strong coupling is shown for representative bosonic and fermionic cases 
in Fig. 4. Specifically, equal-degeneracy pairing models {rii — 712 = ^) are considered, 
and the {V1V2) — (00) states are shown. Here a sufficiently large total occupancy 
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Figure 4. Eigenvalues of tlie bosonic two-level pairing model, with level degeneracies 
ni = n2 = 5 (at left), and fermionic two- level pairing model, with level degeneracies 
ni = n2 = 50 (at right), for repulsive (at top) and attractive (at bottom) pairing 
interactions, shown for the {V1V2) = (00) subspacc, as functions of the control 
parameter ^ between the weak-coupling and strong-coupling limits. All calculations 
are for TV = 50, thus for fl <^ N in the bosonic case and Q, = N (half filling) in the 
fermionic case. The alternative regime in which the bosonic system also has ^ N 
is shown (specifically, for m = 712 = 51 and N = 50) in the inset to panel (a). The 
Hamiltonians H± of (69) and (70) are used in the calculations. 

(^N — 50) is chosen such that the precursors of the phase transitional singularities 
are readily apparent. Spectra are shown for both bosons [Fig. 4 (left)] and fermions 
[Fig. 4 (right)], with repulsive [Fig. 4 (top)] and attractive [Fig. 4 (bottom)] interactions. 
Qualitatively similar spectra in the bosonic and fermionic cases are obtained when level 
degeneracies for the bosonic calculation (fl = 5) arc much less than the occupancy, 
while the degeneracies for the fermionic calculation (Q = 50) are such as to give half 
filling. Then the "envelope" of the spectrum (the range of eigenvalues at a given value of 
the Hamiltonian parameter ^) is essentially identical for the bosonic case with repulsive 
pairing [Fig. 4 (a)] and the fermionic case with attractive pairing [Fig. 4 (d)], i.e., the 
interactions signs which yield a ground state quantum phase transition. Features to 
observe include the essentially constant ground state energy for ^ < 1/5 and downturn 
[from to —1] for ^ > 1/5, an approximately linear evolution of the highest eigenvalue 
from +1 to 0, and a compression of the level density at £^ for ^ > 1/5, a characteristic 
of the excited state quantum phase transition [14]. 

The structure of the eigenvalue spectrum is likewise similar when one compares the 
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hosonic case with attractive pairing [Fig. 4 (c)] and the fermionic case with repulsive 
pairing [Fig. 4 (b)]. This should hardly be surprising. Indeed, when rii = n,2, the 
eigenvalue spectra for Hamiltonians for opposite pairing signs [e.g., Fig. 4 (a) and 
Fig. 4 (c), or Fig. 4 (b) and Fig. 4 (d)] may be obtained from each other, by negation 
of the Hamiltonian and interchange of the level labels 1 and 2, to within addition of 
a c- number function of ^. Therefore, in the present example, the resemblance between 
Fig. 4 (c) and Fig. 4 (b) is a necessary consequence of the resemblance between Fig. 4 (a) 
and Fig. 4 (d). 

The emergence of finite-size precursors to the infinite-A^ singularities associated 
with the quantum phase transition depends not only on but also on the level 
degeneracies rii and n2. An important distinction therefore arises between bosonic 
and fermionic models [14] . For fermionic systems, the total occupancy A^ is limited to 
ni + n2. Therefore, the limit of large A^ can only be taken if the level degeneracies are 
simultaneously increased. Since at full filling (A^ = ni + 712) the spectrum, like that for 
zero filling, is trivial, it is more informative to take the limit N ^ 00 ai or near half 
filling [A^ = + 712) = However, no such restriction arises for bosonic systems, 
and N ^ 00 can be obtained even for fixed level degeneracies. 

Indeed, for the bosonic two-level pairing models, we find numerically that the onset 
of critical phenomena requires N ^ n, not N ^ Q. The evolution of eigenvalues for 
the bosonic system with the same occupancy (A^ = 50) as in Fig. 4(a), but with level 
degeneracies comparable to the occupation ni = n2{= fi) = 49 ~ 50, analogous to 
"half filling", is shown for comparison in Fig. 4 (inset). The eigenvalue spectrum is 
qualitatively different, as compared to Fig. 4(a) or (d), with respect to each of the 
properties noted above, e.g., the ground state eigenvalue is not recognizably constant 
for ^ < 1/5, there is no apparent change in curvature at ^ = 1/5, and closer inspection 
reveals no level spacing compression of the type associated with the excited state 
quantum phase transition. This is already anticipated from the different eigenvalue 
range (~ 3A^^) in the strong coupling limit, obtained in Sec. 5.5. 

Similar distinctions between the large- A^ limit taken with N ^ Q or N Q are 
obtained for the critical scaling properties, which we defer to a more comprehensive 
study. For now, we restrict attention to the basic energy spectra obtained with the 
present transitional Hamiltonian for the general two-level pairing model, and note that 
the spectrum for finite N depends strongly not just on the total degeneracy n — 711 + 712 
of the two levels but on the equality or degree of inequality of the two level degeneracies 
7ii and 712- 

The transitional spectra for two different bosonic models with total degeneracy 
n = 6, and taken with A^ = 10 {i.e., occupation substantially greater than the 
degeneracy), are compared in Fig. 5 (top): the s-b model [ui = 1 and 77,2 = 5) [Fig. 5(a)] 
and the choice of two levels with equal degeneracies {ui = 3 and n2 = 3) [Fig. 5(b) ff]. 
The transitional spectra for the (^1^2) = (00) subspaces in Figs. 5(a) and (b) are similar 

tfFig. 5(b) also corrects a labeling error in the legend of Fig. 6(c) of Ref. [14]. 
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Figure 5. Eigenvalues for the bosonic two-level pairing model, with level degeneracies 
(a) ni = 1 and n2 = 5 or (b) ni = 3 and n2 = 3, and for the fermionic two-level pairing 
model, with level degeneracies (c) ni = 2 and n2 = 18 or (d) ni = 10 and n2 = 10, as 
functions of the control parameter ^ between the weak-coupling and strong-coupling 
limits. All calculations are for N = 10. Eigenvalues are shown only for the lowest- 
seniority subspaces {V1V2), specifically, those with vi + V2 < 2. The Hamiltonian 
of (69) is used for the bosonic calculations and H_ of (70) for the fermionic calculations. 



to each other. Although only irreps of type (OV2) or {IV2) are obtained in the former 
case, more general irreps {V1V2) are possible in the latter case, naturally leading to a 
more complicated spectrum. In particular, it should be noted that the lowest state from 
each subspacc of the form (fiO) approximately tracks the lowest (00) state in energy, 
and that these states are in fact lower in energy than the (00) state everywhere between 
the dynamical symmetry limits [see the lowest curve for (^1^2) = (20) in Fig. 5(b)]. 
It is perhaps not surprising that, given a repulsive pairing interaction, the energy may 
be lowered by breaking pairs within the lower single-particle energy {i.e., increasing 
Vi). In contrast, increasing V2 also enforces nonzero occupation (A^2 > 1^2) of the higher 
single-particle energy level and is therefore not as energetically prefered. 

For the fermionic system, the difference between the transitional spectra for 
near-equal versus highly-imbalanced degeneracies for the two levels is marked. The 
transitional spectra for two different fermionic models with total degeneracy n = 20, 
again taken with = 10 (which now represents half filling), are compared in 
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Fig. 5 (bottom): for the most extremely imbalanced possible choice of degeneracies 
(ni = 2 and n2 = 18) [Fig. 5(c)] and with equal degeneracies (ni = 10 and n2 — 10) 
[Fig. 5(d)]. The quantum phase transition which occurs for equal degeneracies is washed 
out in the limit of imbalanced degeneracies, as is evident in the simple, near-linear 
evolution of the ground state energy across Fig. 5(c). Such an effect may be expected on 
the basis of the Pauli principle. The lower level, of degeneracy rii = 2, easily saturates 
at full occupancy, so that the dynamics are effectively those of a one-level system of 
degeneracy n2 — 18, which does not support critical phenomena as a function of pairing 
interaction strength. 

7. Conclusion 

Although the existence of duality relations between the number-conserving unitary 
and number-nonconserving quasispin algebras for the two-level system with pairing 
interactions is well known, and indeed these relations have proven useful in practical 
calculations for specific special cases of the two-level pairing model, here we have 
sought to estabhsh a systematic treatment of the duahty relations, both for bosonic 
and fermionic two-level pairing models and for arbitrary choice of level degeneracies. 
A principal goal has been to clarify the relationships between the disparate forms 
of the Hamiltonian encountered in the study and application of these models. The 
results arc intended to provide a foundation for a more comprehensive investigation of 
quantum phase transitions in two-level pairing models — including the dependence of 
scaling properties on the bosonic or fermionic nature of the system and on the level 
degeneracies — beyond the special cases conventionally considered, namely, bosonic s-b 
models and fermionic models with equal degeneracy. The duality between orthogonal 
or symplectic algebras and the quasispin algebras is also relevant to the analysis of the 
classical dynamics of the system, through the associated coset spaces [41]. The dual 
algebras yield complementary descriptions involving classical coordinate spaces with 
different dimensionalities [42]. Finally, although the present derivations were given for 
the case of two-level models, they may readily be extended to the SO(ni + n2 + ■ ■ ■) or 
Sp(ni + n2 + ■ ■ ■) algebras associated with multi-level systems, with generators directly 
generalizing those of (14) and (15), for which the quantum phase transitions have been 
much less completely studied. Physical realizations of interest in this more general case 
include the nuclear shell model and descriptions of superconductivity in metallic grains. 
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Appendix. Spherical tensor commutation relations 

When working with angular momentum coupled products of spherical tensor operators, 
it is convenient to consider the coupled commutator [57,72], itself a spherical tensor 
operator, with components given by 

= 5^(aa6/3|c7)[AW,<^]- (A-l) 

a/3 

To clearly set out the identities used in establishing the commutators of the generators 
in Tables 4 and 5, the basic definitions and properties are summarized in this appendix. 
The use of coupled commutation results bypasses the tedious process of uncoupling the 
operators, i.e., introducing multiple sums over products of Clebsch-Gordan coefficients, 
taking commutators of the spherical tensor components, and then recoupling. 

If both bosonic and fermionic operators are to be considered, a consistent set of 
definitions is obtained if the quantity in brackets is taken to be the graded commutator, 
that is, either the commutator or the anticommutator according to the bosonic or 
fermionic nature of the operators. Specifically, 

[AL^\bP] = A^:^Bf - 9,,,BfA^:\ (A.2) 

where 9ab = + if either A or i? is a bosonic operator, and 9ab = — if both A and 
B are fermionic operators. (For the sake of these definitions, it is assumed that a 
bosonic operator has integer angular momentum and a fermionic operator has half- 
integer angular momentum.) The coupled commutator can be written directly in terms 
of coupled products as 

and obeys the symmetry or antisymmetry relation 

The uncoupled commutators of the spherical tensor components may be recovered from 
the coupled commutators, if needed, by inverting (A.l) to give 

[^'a\Bf] = J](aa6/3|C7)[AW, B%^\ (A.5) 

The product rule for coupled commutators is [57, (6)] 

[{A X BY'\cf^ = e d / }(^ X [5,C](^))('^) 

+^dc/}([AC](^)xS)('^)]. (A.6) 

A second application of this identity yields the double product rule needed for evaluating 
commutators of one-body or pair operators, 

[{A X bY'\{C X 

= ^e/lA;[(-)'+'+'+'(-)"+'+'^''{^ e {}{ c ^ X [B,Cf^Y^^ x 

hk 



Dual algebraic structures for the two-level pairing model 



32 



\e g 

+eec0M{-r^'^^^\-f^'^'^'{l ti}{ldl}[Cx {[A,Df'> X E)^'^)]^^)]. (A.7) 

If operators and are creation operators, obeying cannonical commutation or 
anticommutation relations, the canonical commutators are represented in coupled form 
by [57, (10)] 

[ASt]('=)=aWco (A.8) 

and [A, S]^'^) = = 0. The coupled commutator of two one-body operators is 

therefore 

-eabebcecai-f^'^'^^efl HiMc^x bY^^6ad, (a.9) 



as needed, e.g., for the commutators of the generators of U(ni -|- 712). 
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